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Abstract. We study a parabolic equation for finding solutions to the optimal 
transport problem on compact Ricmannian manifolds with general cost func- 
tions. We show that if the cost satisfies the strong MTW condition and the 
stay-away singularity property, then the solution to the parabolic flow with 
any appropriate initial condition exists for all time and it converges exponen- 
tially to the solution to the optimal transportation problem. Such results hold 
in particular, on the sphere for the distance squared cost of the round metric 
and for the far-field reflector antenna cost, among others. 



In this paper wc study a fully nonlinear parabolic flow toward optimal trans- 
port maps between smooth densities on compact manifolds. Let M, M be two 
n-dimensional closed (i.e. compact and without boundary) Ricmannian manifolds 
equipped with the transportation cost c : M x M 9 (x, x) h-> c(x, i)eRU {+00}. 
For example, M = M and c(x,x) = dist(:r, cc) 2 where dist denotes the Riemannian 
distance function of the given metric. As we see the presence of cut-locus for dist, 
the transportation cost may not be smooth, and we denote sing(c) C M x M to 
be the set of points where c fails to be C°° . Let p. p be two smooth probability 
measures on M, M, respectively. (We use p{x),p{x) to express the corresponding 
densities in local coordinates, i.e. p = p{x)dx, p = p(x)dx.) 

One seeks an optimal map T : M — » M, which minimizes the cost functional 



among all maps F which satisfy F#p = p, i.e., J M f{F{x))p{x)dx = f(x)p(x)dx 
for all / G C°°(M). 

It is well known (c.f. [Brl iGMl IMcl IMTWj ) that such an optimal map uniquely 
exists under some appropriate conditions on c,p, p, and it always is associated with 
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a potential function u : M > R, such that the optimal map T (if smooth) satisfies 

(1.1) -D x c(x,T(x)) = D x u{x). 

Here, D x denotes the derivative in the first variable x S M. Note that the condition 
T#p = p forces the potential u (if smooth) to satisfy a Monge- Ampere type equation 

(1.2) det( Uij {x) + Cij (x,T(x)) = \detc fj {x,T(x))\ 

where the subindices i, j, ■ ■ ■ denote the derivative (in local coordinates) in the first 
variable x while the barred subindices i, j , • ■ • denote the one for the second variable 
x. 

The seminal work of Ma, Trudinger and Wang jMTWllTWT] (which follows that 
of Delanoe [D], Caffarclli [Call ICa2) and Urbas [U]) studies regularity of optimal 
maps for general cost functions c by considering regularity of solutions to the fully 
nonlinear (degenerate) elliptic equation (|1.2[) . Importantly, they have identified a 
structure condition on c, now widely called the MTW condition (see 12.11 for more 
details), which later is shown by Loeper [Lolj to be a necessary condition for the 
regularity of solutions to p. 21) . This then is followed by many works including 
those of (in alphabetical order) Delanoe, Figalli, Ge, Kim, Liu, Loeper, McCann, 
Rifford, Trudinger, Villani, Wang and Warren [LaD Eo2 lTW2l iKMll iKrm lKM2l 
rDUl El El EH ILTWl IFKM1I IKMWI IFKM21 IFKM31 IFW] , among others, which 
study both regularity of optimal maps and geometric issues related to the MTW 
condition. 

With this in mind, in the present paper we study the following fully nonlinear 
parabolic equation for u : M x [0, oo) — ^ R, 

Oti 

(1.3) — = lndet(u.y + Cij(x, T(x))) — In p(x) + In p(T(x)) — lndet \cq{x, T(x))\. 

Here, to define T uniquely by u through p. II) and to make the logarithm make sense, 
we assume (as in [MTW]) throughout the present paper, the following conditions 
for the cost function: 

(1.4) the map x i— >• D x c(x, x) £ T*M is injective for (x, x) € M x M \ sing(c); 

(1.5) det Cfj(x, x) ^ on M x M \ sing(c). 

Notice that the right hand side in (| 1 . 3[) is a coordinate invariant quantity. 
The main result of the present paper is the following 

Theorem 1.1 (Parabolic flow toward optimal transport on manifolds). Let M, M 

be two n- dimensional compact Riemannian manifolds without boundary, equipped 
with a cost function c : M x M — > RU {+00} satisfying (|1.4[) and (|1.5[) . Assume 
that c is locally semi-concave on the domain where its value is finite. Let p, p be 
smooth, positive, probability densities on M , M , respectively. Assume further c 
satisfies (i) MTW(<5) condition for some 8 > 0, and (ii) stay- away- from- singularity 
property (see Sections \2.1\ and \2.2\ for definitions) . Let uq : M — > R be a C 2 locally 
strictly c-convex function: see (|2.5p for definition. Then, there exists a unique 
smooth solution u : M X [0, 00) — > R, of (|1.3p with it(-,0) = Uq. Moreover, the 
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derivatives of u (i.e. ||tt||c m for each m £ NJ are bounded uniformly in time, and as 
t — > oo, u(-,t) converges exponentially to a solution of the equation ()1.2j) . and this 
solution to (|1.2[) defines the unique solution to the optimal transportation problem. 

Remark 1. In fact, as a corollary to Theorem 11.11 one can show that C 2 locally 
strictly c-convex functions are indeed globally c-convex in the situation of Theo- 
rem [TTTJ See Corollary [73] (see Section l2~2l for definitions). 

The main consequences of Theorem 11.11 are in the following 

Corollary 1.2. The same existence and exponential convergence result for the 
solution to (|1.3|) as in Theorem ] 1 . II holds for the following cases of M = M, where 
the condition MTW(<5), 5 > 0, and the stay-away-from-singularity property as well 
as regularity of solutions to the elliptic equation (| 1 .2[) are shown in the papers cited 
correspondingly : 

(1) c = dist 2 /2 on the round sphere (see |DL[ ILo2j for elliptic case); 

(2) c = dist 2 jl on small C 4 perturbations of the round sphere, but in a restricted 
sense that Ai,A2 in (|2.4p are restricted by the size of perturbation (see 
[DGl IFRVlj for elliptic case); 

(3) c = dist 2 / 2 on the covering quotient S n /T, T discrete group, e.g. RP". 
(see [DG| IKM2] for elliptic case), and their C 4 perturbations (see [LV] for 
elliptic case). 

(4) c = dist 2 /2 on the complex project space CP™ and the quaternionic n-space 
HP™ with the metric induced from the round sphere by the Hopf fibrations 
S 2n+i _^ cp n md S 4n+3 H pn ^ ee for elliptic case). 

(5) the far-field reflector antenna cost c(x,x) = — \og\x — x\ on the imbedded 
round sphere S n C R™ +1 (Schniirer |Sch| gave a proof in parabolic case, 
see [Lo2j for elliptic case.) 

Remark 2. In all of the above cases one can take uq = as an initial condition. 

Remark 3. Recently, in |FKM3] regularity of solutions to the elliptic equation (|1.2[) 
is shown for c = dist 2 /2 on the multiple products x • ■ ■ x S 1 ™^ of round spheres 
of arbitrary dimension nt and size r,, by showing the stay-away-from-singularity 
property and using the results in [FKMli fLTWj . However, we do not have the 
corresponding parabolic version yet. 

The proof of Theorem 11.11 is given in the rest of the paper; especially, see Sec- 
tions [4] and [Jj We use the tensor maximum principle method for obtaining the 
second derivative estimates, which is essentially the calculation of Ma, Trudingcr 
and Wang |MTWj . given the stay away assumption. We use the parabolic Krylov- 
Safonov theory to obtain C 2,a estimates from C 1,1 estimates, although we note 
that a long and straightforward parabolic adaptation of the Calabi computation in 
[CNSj would furnish the C 1,1 to C 2 ' 1 estimates, if one were interested in a "bare- 
hands" proof. Higher estimates follow by the Schauder theory. In order to obtain 
the exponential convergence, instead of relying on the quite general Krylov-Safonov 
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Harnack estimates, we directly prove Li-Yau type Harnack estimates. Also, a sim- 
ple topological argument is used to show that the resulting limit function at t = oo 
is indeed the solution to the optimal transport problem. 

The proof of Li-Yau type estimates illustrates the natural relevance of the Rie- 
mannian submanifold geometry that is first seen in |KMlj , further refined in [KMWj 
and then again in |Wrj . In this paper, we see that the quantity 8, defined to be the 
right-hand side of (|1.3|) . which measures how far the map at a given time is from 
being volume preserving, satisfies a heat-like equation, which makes it vulnerable 
to the Li-Yau approach. For nonlinear equations, the idea of using the linearized 
operator to define a Laplacian with respect to a metric goes back to Calabi |Cbj . 
and also arises naturally when studying minimal submanifold equations. In our 
case, there is an interesting conformal factor, which forces us to perform a slight 
workaround in two dimensions. See Section [SJ 

Having practical applications in mind, the exponential convergence in Theo- 
rem [LI] can be of independent interest since one can view our parabolic flow as an 
algorithm to construct the solution to optimal transport problem. Indeed, Schnurcr 
[Sch] applied a parabolic flow to construct reflector antennas for given light sources: 
this practical problem in geometric optics has been known to be an optimal trans- 
port problem as found by X.-J. Wang |Wg| . Recently, Kitagawa |Ktj has a result 
similar to ours, which deals with oblique boundary value problem on domains in 
R™, for cost functions satisfying the MTW(O) condition (see Section l2~Tj) . 

The paper is organized as follows. First, after a few preliminaries in Section [2] 
we obtain in Section[3]the stay-away-from-singularity property along the flow which 
is then used together with the maximum principle argument of Ma, Trudinger, and 
Wang |MTWj applied to the parabolic setting, to get uniform second derivative 
estimates along the flow. In Section[4]we use these estimates to obtain the long time 
existence and uniform derivative bounds of the solution. Following some remarks 
about the linearized operator in Section [SJ we prove the Harnack inequality in 
Section [6] Finally, in Section [3 the exponential convergence to the stationary 
solution at infinity is obtained, and this stationary solution is shown to be the 
solution to the optimal transport problem, thus finishing the proof of Theorem ll.fi 
As a corollary to Theorem ll.il it is shown in Corollary 17. II that C 2 locally c-convex 
functions arc globally c-convex. 

Acknowledgement: Y.-H.K. is pleased to thank Albert Chau for helpful dis- 
cussions. M.W. is pleased to thank Jun Kitagawa and Oliver Schniirer for helpful 
conversations. 

2. Preliminaries 

2.1. Ma- Trudinger- Wang curvature. We first explain the MTW curvature that 
is first introduced by Ma, Trudinger and Wang |MTWj as a quanitity which can 
be used to gaurantee interior regularity for solutions of the elliptic equation (|1.2[) . 
This notion is further investigated by Loeper [Lolj and then by Kim and McCann 
[KMT] (see also |Km| lKM2l E3 El EE EE])- As shown by [KMlj . the MTW 
tensor associated to c : M x M —5- R, can be understood as the curvature of a 
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pseudo-metric 

(2.1) h = - (dx l ® dx j + dx j <g> dx^j 

defined on M x M\ sing(c). For (V, V, W, W) G T (x , x) {M x M)x T (x ^(M x M) 
(i.e. V, W G T^M, V,W € T X M), it is computed in local coordinates [MTW] as 

MTW(y, v, w, w) = ( c ifkl ■ r^^'c^n-w 

At this point we mention the (2, 2) form of the tensor, which appeared first in 
[T] . While the coordinate invariance implicit in [T] is easily checked directly by an 
elementary calculation (without refering to curvature), we mention how it comes 
from (|2.ip . Due to the structure of the metric (|2.1[) the sharp operator, which 
identifies covectors with vectors, on M xMisa actually a map from T*M to TM. 
Thus we have the following (2, 2) tensor: 

MTW : TM x T*M x TM x T*M -» K 

(2.2) MTW(F, 77, W, C) = (-Ci jf p + c zjS c Sm c mP r) ^P^O- 

This formulation of the tensor arises naturally in the MTW calcuation and has the 
advantage of having quantities all in terms of vectors and covectors on M. We say 
that c satisfies MTW(<5) condition with respect to a metric g on M if, for S > 

(2.3) MT\N(V, V ,V, V )>S\\Vf g \\ V \\ 2 g 

for all vector covector pairs with rj(V) = 0. The metric g on the source needs nothing 
to do with the cost function, but by fixing the metric, (|2.3p becomes invariant. 

2.2. Stay away from singularity property for a cost function and c-convexity 
of potential functions. Fix Riemannian metrics g, g on M, M, respectively. We 
say that a cost c : M x M — > Ru{+oo} has the stay- away- from- singularity "property 
if for each < Ai, A2, there exists e > depending only on Ai, A2 and c such that 

(2.4) Ai < I det DT\ < A 2 =^> dist(graph(T), sing(c)) > e 

for any differcntiable map T : M — > M given by ([l.lj) with C 2 locally strictly c- 
convex potential function u : M — > R. Here, the C 2 potential function u is called 
locally strictly c- convex if 

(2.5) u ij( x ) + Cjj(x,T(x)) is positive definite for each x G M. 

This local strict c-convcxity on u is equivalent to that the map T is a local diffco- 
morphism. In the above, det DT is computed with respect to the metrics g and 
<7, and the distance function dist is with respect to the product metric g © g on 
M x M. 

To be a solution to the optimal transportation problem for p and p, a solution 
u : AI — > R of (|1.2[) has to be a global c-convex function. Namely, u is given as a 
pair (u, u) as 

u(x) = sup — c(x, x) —u{x), u{x) = sup — c(x, x) — u(x) 
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for all (x, x) £ M x M . If M is a closed manifold and the cost function c is locally 
semi-concave on the set where its value is finite, then for any C 2 locally strictly 
c-convex function u on M, the global c- convexity is implied if the corresponding 
map T via (jTTTJ) is a global diffeomorphism. To see this, suppose u is not globally 
c-convex. Then, there exists xq £ M and x\ 7^ xq such that 

u(xi) + c(x 1 ,T(x )) < u(x ) + c(x , T(x )) 

Thus on the closed manifold M, there is an absolute minimum point zq ^= xq, 
of the function u(-) + c(-,T(xo)). Near zq the cost c(-,T(xq) should be finite by 
the minimum property, thus locally semi-concave by the assumption on c. Since 
u is C 2 , the function u(-) + c(-,y) is semi-concave too, hence supcrdiffercntiablc. 
Thus, this function cannot achieve a minimum at points of nondifferentiability, so 
we conclude that u(-) + c(-,T(xq)) is differentiable at z . Therefore, 

Du(zq) + c(z ,y) = 
Du(x ) + c{x ,y) = 0. 

In particular the local diffeomorphism T is not one-to-one, showing the claimed 
equivalence. 

The stay-away-from-singularity property is shown for the round sphere M = 
M = S n with c = dist 2 /2 by Delanoe and Loeper |DLj . and later also for the 
reflector antenna cost c = — log \x — x\ on S n C R n+1 |Lo2j , and for c = dist 2 /2 on 
the perturbations of the round sphere and its discrete quotient jDGl ILV1 IKM2) , on 
the Hopf fibration quotients of the sphere such as CP" and HP" KM2 . Recently, 
such result is also shown for c = dist 2 /2 on the products of round spheres M — 
M = S% x • • • x of arbitrary dimension and size |FKM3j . 

Remark 4. Strictly speaking, the stay-away results in the papers |DL[ ILo21 IDGi 
IhVl IKM2[ IFKM3] are shown with respect to globally c-convex functions, which in 
general differ from locally c-convex ones. However, in these cases, one can actually 
prove such stay-away property with respect to locally strict c-convex functions. 
Alternatively, it can be shown that in those cases, locally c-convex functions are 
globally c-convex, using the results of (TW2l El ITW31IFRV] . 

3. Estimates 

In the following, T will always denote the map given by the relation (|l.ip . which 
is equivalent in any coordinate chart x — (x , ■ ■ ■ ,x n ), to 

(3.1) u t {x) + Cl (x,T{x)) = 0. 

Further differentiation gives 

Uij + Cij(x, T) + Cis{x, T)TJ = 0. 

Here and henceforth we use f,s,p,v, etc to denote differentiation in the second 
variable of c(x, x), with the map T represented in coordinates as T s , s = 1, • • • , n. 
Taking the determinant of the above equation gives the elliptic equation (|1.2j) . A 
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useful observation is that the matrix Wij(x) := + Cij(x,T(x)) (when positive 

definite) gives a Riemannian metric on M by the following identity: 

where (Id x T)*h is the pull-back of h by Id x T : M 3 x i-> (x, T(x)) G M X M. 
For future reference we note the following results of differentiating 



(3.2) Tj 



(3.3) = -(#— Wij + c l fc ls (c lpk + ci w Tl) Wi j . 

Here, to be clear, we write -^-Wij to express coordinate derivatives of this tensor. 
Subscripts for other functions mean the corresponding coordinate derivatives. Su- 
perscripts such as c ls denotes the (i, s)-entry of the inverse matrix of ^^ a _ t cj , etc. 
Defining 

(3.4) 6(x,u) = lndet my - \np(x) + lnp(T(x)) - In | det c^x, T(x))\, 
the flow (|1.3[) is rewritten as 

(3.5) ut = e. 

Here and henceforth, the subscript t denotes the time derivative. 

3.1. Stay away from singularity. We make use of the linearized operator L of 
at a function u. This operator L is covariant, and as seen from [TW] the operator 
L has the following local expression. 

(3.6) Lv = w ij Vij - (/r'"V,.,, r - ,v - ;,,) 1 /',«'* r - c iS c m c pk ) v k . 

where all coefficients are computed at (x, T(x)). It is important to notice that there 
is no zeroth-ordcr term in L so that we can apply parabolic maximum principle. 



Lemma 3.1. If u is a solution of \3.5\) . then 6 satisfies 
(3.7) 6 t = L6. 

Proof. This follows from the definition of linearized operator (8(x,u)) = L-j^u 
and using the equation (|3.5p . □ 

Proposition 3.2 (Stay-away-from-singularity property along the flow). Let M, M 

be two n-dimensional compact Riemannian manifolds without boundary, equipped 
with cost function c : M x M — > R satisfying (jl.4|) and (jl.5p . Let p, p be smooth, 
positive, probability densities on M , M , respectively. Suppose the cost c satisfies 
the stay-away-from-singularity property (see Section ] 2. There exists a constant 

e = e(e| t=0 ,l|lnp|| OO) ||lnpl| oo ) > 0, 

such that in the time interval of existence of the smooth solution u to the para- 
bolic equation (|1.3j) . dist(graph(T) , sing(c)) > e for T given by (|l.lj) . Here, £ is 
independent of the time t. 
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Proof. It follows from Lemma 3.1 and the parabolic maximum principle applied to 
(|3.7p that 9 is bounded along the flow. Now as 

6(x,t) = lndetDT - ]np(x) - ]np(T(x)), 

a bound on || lndetDTHoo follows from the bound on 6. This together with the 
stay-away-property (|2.4p completes the proof. □ 

Corollary 3.3 (Uniform bounds on derivatives of c along the flow). Suppose the 
same assumption as in Proposition \3. 2\ holds. Then, in the time interval of existence 
of the smooth solution u to the flow of (|1.3j) , each derivative of c in any order such 
as Cij, Cis, Cijki c ijki etc. a t { x i T(x)) computed in any fixed coordinates, is uniformly 
bounded, especially independent oft. Moreover, Ci S is uniformly away from 0. 

3.2. Second derivative estimates. On the both source and target manifolds, we 
work in normal coordinates throughout this section. Noting that the transition 
functions between charts change in a predictible way, an e/4 covering argument 
can be used to show there is a uniform bound on (in particular third and fourth) 
derivatives of the cost function which doesn't depend on the charts, regardless of 
the point at the origin. Thus throughout this subsection, we assume the result of 
Corollary 13.31 which is an essential ingredient in the following estimates. Following 
the calculations given in |MTWj we show 

Lemma 3.4. Assume the uniform bounds in Corollary \3.S\ Let u be a solution of 
h3. 5\) . and let v be a coordinate direction in a local chart. Then 

d 

(3.8) ~Qi Wvv ~ ( Lw ^ vv - ~ w%3 ( <"'.//" + c ijs c sm c m f P ) c pk c rl W[ v w k v 

+ C{p, p, D\ g)(l + E w u + E w " E 
Here L is operating on the (2, 0) tensor using metric covariant differentiation. 

Proof. To simplify the calculation choose normal coordinates around x such that 
v = gfp Now from (|3~1J) . 

(Lw) u = w ij {wxi, ij -w ij>11 )+w ij w ii>11 -{w ij c^sC^ - {p^pfC^ - c"c Mp c pk ) d k wu, 
with 

mj,u = Uijn + Cijn + 2cijsiTl + CijsT^ + c ilW TlT{ + w * d 2 g 

wu,ij = uiuj + ciuj + cumTj + cusjT? + cusT^ + cupfTfT- +w* d 2 g 

where w * d 2 g denotes a linear combination of second derivatives of the metric with 
components of w. Thus 

(Lw) n > -K + w l3 {c lls T° 3 + c llw T?Tr - c, ,j ; : - c^TfTf) 

+ w^dnWij - (w tJ c ijs c sk - (\np) f c kf - c is c m c pk ) d k w xl . 

Here as in |MTW] we let 

K = C E w " E w n + ^E' r " • r E"- 2 <- 
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Now use ()3.3[) to see 

(3.9) T n = -c ks w khl + c k fc ls (c lpl + c m T[))w kl 
and 

(3.10) w^Tlj = -c ks w^ {d jWki - d kWij ) - 'a-",),,,;, + c^c ls {c lpj + c tpf TV). 
Note that 

djWk% - dkWij = CkisT- - c ijs T k 
and use this to cancel one dkWn term to get 

(Lw) n > -K - cn s c ks w l ' J dkWij+w lJ ' (-Cij S cipfT[wki - Ci jP rT^T[) 

+ w ij d n w l3 + ((lnp) f c feF + c lS c m c pk ) d k w u . 

Now using 

n = w^dnwii + dyUhwij + Qnp)^ + c^c^Tf, + r'>,.. ; „ 7f/ - K 
> w ij d 11 w ij + ((lnp)-+c is c ig p) c ks w kltl - K 
we have 

Mn > -K - c lls c ki 6 k + flu 

w ^ ^ CijsC cipf -\- Cijpf^j T^T^ 

Now, finally we compute 

d 

— wn = 6»n + ciijT/. 
at 

Differentiating (|3.1[) with respect to t we get 

u k t + c ks T t s = 0, 

which, subtracting, we arrive at 

(Lw) n - wnt > -w %3 (~c ljs c ls ci P f + c ljpf ) T[T\ - K 

which is the conclusion of the Lemma using (|3.2p . □ 

Theorem 3.5. For a closed Riemannian manifold M and an interval [0, 1), suppose 
that u : M x [0, 1) — > R is a smooth solution to the parabolic equation (|3.5p . Assume 
that c satisfies (i) the MTW(<5) condition, and (ii) the stay- away- from- singularity 
property (so that the uniform bounds of Corollary Iff. 51 hold on M x [0,1)). Let 
v G TM, \v\ = 1, be such that 

w vv = max w zz . 

z£TM,\z\=l 

There is a constant C = C(6,c,n,g), especially independent oft, such that if 
w vv {x, t) > C then 

w vvt {x,t) < 0. 

In particular, w zz {x,t) < C for any (x,t) £ M x [0,1) and \z\ = 1. 
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Proof. We use Hamilton's jHaj parabolic maximum principle argument for tensors. 
Analyze the first term in the right hand side of (|3.8j) 



(3.11) -w^MT\N%wi„w kv 
where 

Diagonalize w with v = and (|3.11|) becomes 

•t>i 

By the MTW (S) condition, we have 



2 
9 

i>l i>l 



thus p. lip is bounded by 

(3.i2) c J2 m - 6 w " w ^ ^ w » - Scn J2 '• 

i>l i 

as our chart for the source M is normal at x$. Since detwy is bounded it follows 
from the arithmetic-geometric mean that 

i/(«-i) 



So finally, from Lemma 13.41 



d_ 

dt Wvv 



{Lw) vv w u - 5c n Y ™ U Wii ■ 

i 

i 

+ CY W "(Y1 W H ~ | c « w ii) 



2 

- _ 2 C " W ") C " (51 Wli ) + C ' 



We see that when ^ wu is sufficiently large, the right-hand side must be negative. 
This completes the proof. □ 

Corollary 3.6. In the situation of Theorem Iff. 51 the spatial second derivatives of 
the solution u to the parabolic equation (| 1 . 3[) remain bounded (uniformly in time) 
and u stays locally strictly c- convex. 

Proof. An upper bound on the eigenvalues of u>y is given in Theorem 13.51 From 
the identity u>y = Uy + cy, the bound on ity follows. Because 9 remains bounded 
by parabolic maximum principe for 9 t = L9, we have a lower bound on determinant 
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det Wij. An upper bound on eigenvalues, plus lower bound on determinant implies 
lower bound on the eigenvalues. It follows that 

for some e > 0, thus local strict c-convexity of u follows. □ 

4. Proof of Theorem 11.11 Existence of solution and uniform bounds 

In this section we show that the solution to parabolic equation (|1.3j) exists for 
all t 6 R + under the assumptions of Theorem 11.11 We also show the solution 
has uniform C k derivatives in x, where each C k norm is uniform in t. Through 
this section we use Corollary 13.31 in an essential way. In the following, all the 
C 1 ' 1 , C 2 ' Q , C 4 ' Q estimates and so on, are estimates on the derivatives in x, and are 
all uniform in the time variable t. But, by C a we will mean C a both in x and t. 

Short-time existence: Since M is a closed manifold, a standard theory implies 
the existence of a short-time solution to (|1.3|) for any locally strictly c-convex smooth 
initial data, and that from Corollarv l3.6l the solution u is locally strictly c-convex on 
the time interval of existence. (See jKtj for a proof of short time existence regarding 
the same equation with a more involved boundary condition.) 

Long-time existence: Apply Thcorcm l3.5l (or Corollarv l3.6[) to get C 1,1 estimates 
for u. This in particular makes the equation (|1.3p as well as the linearized equation 
v t = Lv (see (|3.6|) ) uniformly parabolic with bounded coefficients. Now applying 
Krylov-Safonov theory (c.f. |Lb| Lemma 14.6]) to (|1.3[) one has C 2 ' a estimates. 
From the short-time existence above and Arzela-Ascoli, this shows that the solution 
cannot cease to exist at a finite time, thus exists for all t £ [0, oo). 

Uniform C k bounds: To sec uniform C k bounds, first differentiate (|1.3[1 with 
respect to any coordinate direction ^pr, to see satisfies 

(4.1) uh = Lu k . 

From C 2 '" estimates of u, this linear equation for u k is uniformly parabolic with C a 
controlled coefficients (here C a control in t follows from the parabolic equations like 
(|4.ip with space C 2,a estimates), and in particular, we can apply parabolic Schauder 
estimates to conclude that u k has C 2 ' a estimates, thus obtaining C 3 ' a estimates for 
u. Similarly, differentiating (|4.1[) we obtain a parabolic equation for u^i also with 
coefficients and inhomogeneous terms all of which arc C a , and thus follows C 4,a 
estimates of u, so forth. Thus, we have uniform C k , k = 0, 1, • • • , bounds for u as 
claimed. 

To obtain the exponential convergence, we need a Harnack inequality, which is 
shown in the next sections. 
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5. A Li-Yau type Harnack inequality 

In this section, as a preliminary step to the proof of the exponential convergence 
to a solution of the elliptic equation (|1.2p , we derive a Harnack inequality for the 
quantity 6. We first find a simple expression of the linearized operator L of 6. This 
expression allows us to derive a Harnack type estimate (see Theorem 15. 2|) . whose 
corollary (Corollarv l5.3j) is used to show the exponential convergence in Section [7] 

Let us find a convenient expression for the linearized operator L, using jWrl Prop 
2.1]. As a linearized operator of coordinate invariant fully nonlinear equation Q1.3p . 
L is expected to be related to a Laplacian operator of certain Riemannian metric. 
More precisely, for a manifold M, suppose that a map T is c-expdu for some u. 
(This means by definition T is given by (|1.1|) .) Define a function tp by 

_ ( p\T(x)) dot DT \ mn ~ 2) 
[ ' W \\detc lS (x,T(x))\) 

We observe the following 

Proposition 5.1. Letn > 3. The linearized operator L of 9 (see (|3.4|) ) is expressed 
as 

Lv = tpA^ w v, 

where A^ w is the Laplace- Beltrami operator with respect to the metric ipWij given 
on M. 



Proof. First of all, let 

/ p{x)- P (T{x)) y /( - 2) 



(5.2) 9ij{x) = Wij(x) 

Also, recall 



\\det Ci S (x,T(x))\ 



6 = lndct DT - In p{x) + In p(T(x)). 
Then [Wrl Prop 2.1] states that 



( ^mx^v'^ ( i 

\detc ls {x,T{x))) {^ 9V+ 2 



A general formula for conformal metrics shows that if g = e'g, then 



~ 77,-2 

s'A^ = A0+-^-(V^,V/>. 



It follows immediately that 
Lv 



( p{x)p{T(x)) V /( "' 2) e/(n _ 2) 
Vdetc lS (x,T(a;))/ 
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Using the expression for 9 we observe that 

( p{x)p{T{x)) \ 1/[n ~ 2) e/(n-2) 



\detc is (x,T(x)) 

( p(x)p(T(x)) \ 1/(n ~ 2) f(detDT)p(T(x))\ 1/(n - 2) 



\detc lS (x,T(x)) J \ p(x) 

p 2 (T(x)) det DT^ 1/( "~ 2) 



det c is (x,T(x)) 

The result follows. □ 



Noting that by Lemma \3.1\ and Proposition [53] 

e t = L9 = i'A^e 

we derive a Harnack estimate for the operator L, when n > 3. Here, the expression 
ip/S.^ w enables us to easily modify the argument in jLYj to obtain 

Theorem 5.2 (Harnack inequality). Let M be a compact manifold of dimension 
n and let g(t), < t < oo be a family of Riemannian metrics on M such that 

(5-3) -J- 5 (0) < g(t) < C o g(0) 

d 



dl 9 



< Co 
%(*)> -K 9ij (0) 



with universal constants Co,K > 0. Let X(x,t) be a positive function with deriva- 
tives uniformly controlled (independent of(t,x)) and bounded above and away from 
zero. Let U(x, t) be a nonnegative solution to 

(5.4) U t = X(x,t)A g(t) U. 

Then there exists a constant C > depending only on Cq,K, <?y(0) and the bounds 
on the derivatives of A so that for < t\ < t 2 < oo 

sup U(x, t x ) < inf U(x, t 2 ) C— exp <^ -2- + C(t 2 - h) 

xEM X ^ M tx I (t 2 - tx) 

Note that by the a priori estimates in the previous sections, the metric gij = ipWij 
and scalar A = ij) all satisfy the assumptions in Theorem 15 .21 thus we obtain the 
corresponding Harnack inequality for the operator L. In particular, we have 

Corollary 5.3. For n > 3, let U : M x [0, oo] — > R be a solution to the parabolic 
equation U t = LU , where L is the linearized operator in (|3.6[) . There exists a 
constant independent of (x,t) £ M x [0,oo] such that 

sup U(x,t + 1/2) < C inf U(x,t). 

ism X ^ M 
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6. Proof of Theorem 15.21 (Harnack inequality) 

This whole section is devoted to the proof of the Harnack inequality claimed in 
Theorem 1 5. 2 1 Since the equation (|5.4p is not a pure heat equation, but conformally 
related to one for a time dependant metric, we are forced to reprove the Harnack 
estimate for this operator. The following argument is a slight modification of that 
found in |LYj . Let 

f = togU, 



then 



Now let 



ft = A(A/ + |V/| 2 ). 



F = t(X\Vf\ 2 -af t 



We directly compute 



AF > t 



AA|V/| 2 + 2(VA,V 

—a 

2 



V/| 2 )+AA|V/| 2 - a (A/) t 
(|V 2 /| + |V/|) 



> ^ / 2A||V 2 /|| +2A(V/,VA/} -2XK\Vf\ 2 -e||V 2 /|| - ±f\Vf 
~ { -C 2 |V/| 2 -a(A/) t - £ ||V 2 /|| 2 -f 

> t (7 || V 2 /f + 2A (V/, VA/) - a (A/) t - C 4 |V/| 2 - C 3 



In the second line we applied the Bochncr formula, with K a uniform lower bound 
on the time-dependent Ricci curvature, C\ is an upper bound on ||VA|| , e is 
somewhat smaller than a lower bound on A, C2 bounds AA and the constant C3 
bounds the time derivative of the metric. The constant 7 is a positive lower bound 
for A. 
Now 

(6.1) A/ = -|V/| 2 + l -U = ~ (jF (1 - a)/ t ) . 



We thus further estimate 



AF> t 7||V 2 /|| 2 -2A(V/,Vl(^-(l-«)/0) 
" V a[$(±F-(l-a)ft)] t -Ci\Vf\ 2 -C 3 



> t7||V 2 /|| -2(Vf,VF}+2t(l-a){Vf,Vft)-C 5 t\Vf\\(XAf)\ 

1 / F 
+ -laF t - aj 



'aF t -a^-- ta(l - a) ftt ) ~ C 6 t\(XAf)\ - C 4 i|Vjf - C 3 t 
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> t l2 || V 2 /f - 2(V/, VF) + 1 {2iA(l - q)(V/, V/*) - to(l - «)/«} 
+ i ^F, - a|) - C 7 t\\7.f\ 2 - C s t 

= tl2 ||V 2 /|| 2 - 2(V/, VF) + 1 1(1 - a)F t - (1 - a) j - AtC 3 |V/| 2 
1 ( „ F\ „ |2 



_ _ a _ j _ C 7 i|V/| 2 - c 3 t 

> *72 ||V 2 /|| 2 - 2(V/, VF) + If* - ~ - iC 9 |V/| 2 - C 3 i 

> (-A|W| 2 + f t ) 2 - 2(V/, VF) + if! - ~ - tC 9 |V/| 2 - C a t. 

Here C5 bounds |Vj|, and Cq bounds |jL, with 72 an even smaller constant, 
related to C7. In the last line we applied || V 2 /|| 2 > ^(A/) 2 and the relation (|6.1|) , 
Now at a maximum for F which happens at some positive time on M x [0,T] 
we can conclude 



> t^- (-A|V/| 2 + /) 2 - 2(V/, VF) + ~F t - ™ - iC 9 |V/| 2 - C 3 t 
An v ' A A t 

(6.2) > t 73 (-A|V/| 2 + / t ) 2 -if -tC 9 \Vf\ 2 -C 3 t. 
Now let 

y = A|v/| 2 , 
« = /*. 

Continuing to follow jLYl Eq 1.9] we expand 
(6.3) 

{y-zf =(--£) ( 2/ -az) 2 + (l-e-5-i + ^)?; 2 +(l-a+^a 2 )z 2 +£y(y-az)+(5y 2 
\ a 2 / a 2 2 



Choose 



s 1 



and a so that 



1 

- > e > 

a 

*>!>o. 

Specifically, one may choose a = | to satisfy these conditions. Note that the second 
and third terms in (|6.3p vanish, and multiplying (|6.2[) by i, absorbing bounds on A 
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and combining with (|6.3 



(I. £ ) ,, +74 (,_•) tV - F -e [Jgffi + 



> 74 (-- £ J F +74^-tt)<V-^-^ 
From which wc conclude that 



-Cut 2 . 



F < 



ClO 



C 74(5-f) 



274 (i-e) 7 4 (i-e) 



(6.4) A|V/| 2 -a/ t < y + C. 

Now we have arrived at this conclusion assuming that the maximum happens away 
from t = 0. Letting / = log({7 + t) this assumption is available. We then take 
i -> 0. 

Now consider now a path 7 : [0, 1] — s- M x [ii,f2] such that 7(0) = {y,t 2 ) and 
7(1) = (x, ti), which for convenience, projects to a geodesic in M tl and has constant 
speed in t. Using the assumption (|5.3[) and (|6.4D 

f(x,h)-f(y,t 2 )< f {\Wf\C d(x,y)-(t 2 -h)(f) t }ds 
Jo 



C 



^ 177 ^—^C 2 d 2 (x,y)+C(t 2 -t 1 ) + Clog^ 

4(t2 — ti) mm A ti 

where in the last line we integrate using that t = (1 — s)t2 + sii. This completes 
the proof of Harnack inequality (Theorem I5.2j) . 



7. Proof of Theorem 11.11 Exponential convergence to the solution 
to optimal transport problem 

This section completes the proof of Theorem 11.11 Wc assume the conditions in 
Theorem 11.11 In the first subsection we show the exponential convergence to a 
stationary solution, and then in the last subsection we show that the stationary 
solution is indeed the solution to the optimal transport problem. 



7.1. Exponential convergence. 
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7.1.1. Case n > 3. With this parabolic Harnack inequality at hand, the claimed 
exponential convergence to the optimal transportation map follows from a rather 
standard argument. To see this, consider for ieN. 

U k = sup 9(x, k - 1) - 9(x, (k - 1) + t), 
L k =9(x. k-l + t)- inf 6(x, k - 1) 

that are also solutions to the equation U t = LU . By the strong maximum principle, 
both U and L are positive functions for positive t, for all k. Also, let H(t) = 
M0(x,t). The Harnack inequality in Theorem [5T2l yields 

sup9(x,k- 1) -inf0(>,fc- -) < C(sup6(x,k- 1) - sup0(x, &)) , 

sup 6(x, k - -) -inf 0(z,fc - 1) < C (inf 9(x, k) -M9{x,k- 1)) 
for some fixed constant C > 1. It follows by adding the equations together that 
H(k - 1) + H(k - -) < C (if (fc - 1) - iJ(fc)) 

which implies 

ff(fc) < eiJ(fc - 1) 
where e = < 1- By induction we observe 

H{k) < e k H(0). 

It follows that H(t) < Ce-P* where e = e _/3 . Therefore, 6» converges to the limit 
6*00 = const, exponentially fast as t — > oo. Now the quantity 9 can be larger than 
somewhere, only if it is smaller than somewhere, as can be seen by integrating 
the change of measures and using J M p(x)dx = p(x)dx. It follows that 9oc = 
and thus from ut = 9, we see ut — > exponentially fast as t — > oo. This implies u 
converges exponentially fast to a stationary solution Uoo . which is smooth from the 
uniform C m , m € N, estimates. Because 9^ = 0, solves the elliptic equation 
(|1.2[) . Considering the discussion in Section [7?2l this solution is a solution to the 
optimal transportation problem. This finishes the proof of the claimed exponential 
convergence and of Theorem II .11 for n > 3. 

7.1.2. Case n = 2. On M 2 xM 2 a solution u will satisfy all of the estimates, hence 
exists and enjoys subsequential convergence at infinity. The only missing piece is 
the Harnack inequality. However, we can fake a third dimension and get a solution 
on M 2 x S 1 — > M 2 x S 1 by letting u(x, z) = u(x) and taking products of the 
measures with uniform measures on S . Hence u(x, z) will also be a solution, and 
will converge in the same way to the three dimensional product solution. 

7.2. The limit stationary solution is the solution to the optimal transport 
problem: strict global c-convexity: To conclude that the limit stationary so- 
lution, say Uqo, is a solution to the optimal transport problem (and not a spurious 
solution to the elliptic equation), it remains to show that u^o is globally strict 
c-convex, which from the discussion in the middle of Section 12.21 follows if the cor- 
responding map, say Too, is a global diffcomorphism. (This map Too is already a 
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local diffcomorphism by local strict c-convexity of Moo-) To see this, we use that 
Moo satisfies the Monge- Ampere type equation (|1.2|) . If is not one-to-one, as 
a local diffcomorphism between closed manifolds, it is a covering map, having the 
topological degree greater than 1. Thus, from (|1.2[) the push-forward Too#p is a 
multiple of the target measure p. But, this contradicts the fact that J M p = Jv> p 
since /^Too#jO = J M P- This finishes the proof of the fact that the limiting sta- 
tionary solution of (|1.3j) is the solution to the optimal transport problem, and thus 
together with all the previous sections (especially Section 2]) it completes the proof 
of Theorem O 

As a final remark, we state a corollary to Theorem in particular this last 
paragraph: 

Corollary 7.1. Assume that the same conditions as in Theorem \l.l\ hold. Then, 
any locally strictly c-convex C 2 function uq is in fact globally c-convex. 

Proof. Let To denote the corresponding map to uq by the formula Local 
strict c-convexity implies that To is local diffcomorphism and for global c-convcxity 
of uq, it is enough to show that To has topological degree 1. From Theorem 11.11 
the map T depends continuously on the time variable t. In particular, the degree 
stays constant. From the result of Theorem 1 1 . 1 1 the map Too hi the lim t — > oo is a 
diffcomorphism and so its degree is 1. This shows that the degree of To is also 1, 
and completes the proof that uq is globally c-convex. □ 
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